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Here we present Spitzer Space Telescope imaging of Cyg A with the Infrared Array 
Camera at 4.5 //in and 8.0 /xm, resulting in the detection of the high-energy tails or 
cut-offs in the synchrotron spectra for all four hotspots of this archetype radio galaxy. 
When combined with the other data collected (and re-analyzed) from the literature, our 
observations allow for detailed modeling of the broad-band (radio-to-X-ray) emission for 
the brightest hotspots A and D. We confirm that the X-ray flux detected previously 
from these features is consistent with the synchrotron self-Compton radiation for the 
magnetic field intensity B ~ 170 //G in hotspot A, and B ~ 270 /iG in hotspot D. We 
also find that the energy density of the emitting electrons is most likely larger by a factor 
of a few than the energy density of the hotspots' magnetic field. We construct energy 
^ ■ spectra of the radiating ultrarelativistic electrons. We find that for both hotspots A 

and D these spectra are consistent with a broken power-law extending from at least 
100 MeV up to ~ 100 GeV, and that the spectral break corresponds almost exactly to 
the proton rest energy of ~ 1 GeV. We argue that the shape of the electron continuum 
most likely reflects two different regimes of the electron acceleration process taking 
place at mildly relativistic shocks, rather than resulting from radiative cooling and/or 
absorption effects. In this picture the protons' inertia defines the critical energy for the 
hotspot electrons above which Fermi-type acceleration processes may play a major role, 
but below which the operating acceleration mechanism has to be of a different type. At 
energies > 100 GeV, the electron spectra cut-off /steepen again, most likely as a result 
of spectral aging due to radiative loss effects. We discuss several implications of the 
presented analysis for the physics of extragalactic jets. 



Subject headings: galaxies: jets — galaxies: individual(Cygnus A) — acceleration of 
particles — radiation mechanism: non-thermal 
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Introduction 



Cyg nus A is one of the best studied examples of Fanaroff-Riley type II radio galaxies ('classical 
doubles' , iFanaroff & Rilevl 19741 ) , being unique ly observed at different scales and wavelengths t hanks 



to its proximit30 and high total radio power ( Carilli Barthel 1996 ; Carilli fc Harris 1996 ). The 



distinctive hotspots near the extremities of the radio struct ure are roughly 70 kpc f rom the nucleus 
in the North- West (A and B on the jet side; designation by lHargrave Ryldll974l ) and the South- 
East (D and E on the counter-jet side) directions, and are the brightest known radio hotspots. They 
have been detected across the entire accessibl e radio band from < 100 MHz up to > 100 GHz (e.g., 
Carilli et al.lll99ll : IWright &: Birkinshawil2004l : lLazio et al.ll2006l ). These features are understood as 
the terminal regions of relativistic jets, where bulk kinetic power transported by the outflows from 
the active center is converted at the strong shock (formed due to t he interaction of the jet with the 
ambi ent gaseous medium) to the internal energy of the jet plasma (jBlandford Reeslll974l ; IScheuer 
1974 ). The nonthermal hotspot emission is believed to originate mostly in the downstream region 
near the reverse shock front, just before the shocked jet plasma is turned aside by the ambient 



medi um and forced to flow into the extended radio lobe (see a general review by iBegelman et al 
19841 ). 



Besides this general picture, details of both hotspots' macroscopic structure and microscopic 
processes taking place therein are not well explored, mostly because crucial parameters of the 
jets (like the jet composition or bulk velocity on large scales) are still hardly known. Numerical 
hydrodynamical studies indicate a complex and time-dependent morphology of the jet termination 
regions, with an imp ortant role played by turbulent vorti c es and small-scale oblique shocks forme d 



at the edges (e.g., Marti et al 



1997 



Scheck et al. 



2002 



Saxton et al. 



2002 



Mizuta et al.1 12004 ). 



Overall hotspots' structures like the ones observed in classical doubles can be however relatively 
well reproduced in these simulations with the jet outflow assumed to be dynamically dominated by 
cold particles (either electron-proton or pure electron-positron pairs). Much worse is the situation 
regarding microscopic processes leading to the conversion of the jet kinetic energy to the internal 
energy of the jet plasma (both particles and the magnetic field) at the reverse shock front, involving 
particle acceleration and possible energy equilibration among different plasma species. In particular, 
although the diffusive (1st order Fermi) shock acceleration was widely assumed in the past as 
the process of choice fo r generation of the broad-band non-thermal electron en ergy distribution 
within the hotspots (see lHeavens fc Meisenheimeiill987l ; lMeisenheimer et al.lll989l ). the most recent 
studies indicate serious problems with such a picture. These are mostly due to the expected 
(at least mildly) relativistic velocities of the jets in FR II sources on large scales and oblique 
magnetic field configuration, prec l uding particles from undergoing lst-order accel e ration process 



20061 : iLemoine et al.ll2006L discussing different aspects of this issue). 



(see, e.g.. IBegelman Jz Kirkl Il990l : iHoshino et al.1 LL992J; iNiemiec h Ostrowskil [2004J; iNiemiec et al 



1 At the redshift of z = 0.0562 ([Stockton et al.ll 1994 ) the luminosity distance to Cyg A is d L = 248.2 Mpc, and the 
conversion scale is l" = 1.079 kpc for the modern cosmology with Qm = 0.27, Qa = 0.73, and Ho = 71 kms" 1 Mpc" 1 . 
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Some of the problems mentioned above regarding understanding of the hotspots (and thus of 
extragalactic jets) can be overcome by putting observational constraints on the theoretical models 
by means of accessing, and detailed modeling of, the broad-band hotspots' spectra. Cyg A is of pri- 
mary importance in this context, since its bright hotspots c an be detected and resolved at different 
wavelengths. And indeed, ROSAT observations reported bv lHarris et al.l (|1994l ) revealed relatively 
intense X-ray emission of hotspots A and D, being consistent with the synchrotron self-Compton 
('SSC') model for the hotspots' magnetic field of ~ (100 — 300) /iG in energy equipartition with 
the emitting electrons. These important findings — in fact the very first observational sugges- 
tion for the minimum energy condition fulfill ed in the case of ex tragalactic jets — were confirmed 
recently by the Chandra X-ray Observatory (jWilson et al.ll2000i ). Still, the insufficiently sampled 
synchrotron continuum of the analyzed features (unknown minimum an d maximum synchrotro n 
frequencies in particular) made these conclusions not fully robust (see, e.g.. lKino &: Takaharall2004l ) . 
Hence an observational effort was made to cover the whole synchrotr o n spectra in Cyg A hotspots, 



both at th eir low-frequency ( Muxlow et al 



frequency (IMeisenheimer et al.1 119971 : iNilsson et al.1 119971 : ICarilli et al 



1988: 



Carilli et al 



1991 



Laz io et al.1 120061 ) and high- 
19991 ) "ends , as well as to 



construct a large sample of other hotspots with multiwavelength data avai lable (see iBrunetti et al 



20031 : lHardcastle et al.ll2004l : iKataoka Stawarzll2005l : ICheung et al,ll2005l . and references therein). 
This is also crucial for understanding particle acceleration taking place at relativistic shocks. 

Here we present Spitzer Space Telescope imaging of Cyg A with IRAC at 4.5 jum and 8.0 /iin, 
resulting in detection of the high-energy tails or cut-offs in the synchrotron spectra for all four 
hotspots of this archetype radio galaxy. This, together with the archival data collected from the 
literature (see § 2), allows us to construct the broad-band spectral energy distributions for the 
discussed features, to model them in terms of synchrotron and SSC processes (§ 3), and thus to put 
meaningful constraints on the physical processes taking place at the terminal shocks of relativistic 
jets (§ 4). 



2. The Data 
2.1. Spitzer Space Telescope Observations 



We observed Cygnus A with the Infrared Array Camera (IRAC; iFazio et al.l 120041 ) aboard 
Spitzer on 21 October 2005, and the resulting image is shown in Figure 1. Individual 30 second 
exposures were obtained simultaneously at 4.5 and 8.0 fim over a 36 position Reuleaux pattern 
(1/4 th sub-pixel sampling) . This dither pattern results in even coverage of a central ~ 4 x 4 arcmin 2 
field which safely encloses the full extent (~ 140") of the Cyg A radio source. The total field 
coverage is ~ 6 x 6 arcmin 2 with the exposure decreasing by about x 1/4 th of the total in the outer 
~ 1 arcmin wide edges. We repeated the dither 12 times for a total integration time of 3.6 hrs in 
the central portion of the field. The 3cr point source detection limits in the central field are 0.4 
(4.5 //m) and 3.0 (8.0 /im) /iJy (Spitzer Observer's Manual). 
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The pipeline processed basic calibrated data (BCD) files of the individual exposures were ob- 
taine d from the Spitzer Science Center and combined with the MOPEX package ([Makovoz &: Marleau 
20051 ) resampling to 1/3 rd of the native pixel sizes' 1.22" (there are < 1% distortions across the 
FOV and between detectors which we ignore). The hotspots were expected to be faint in the 
mid-IR so we required precise image registration with the other wavelength data (§ 2.2, § 2.3). To 
facilitate this registration, the frametime of 30 seconds was chosen to be long enough to minimize 
the noise per exposure while also short enough to ensure an unsaturated image of Cyg A in at least 
one IRAC band (4.5 /mi). Based on the radio core and position of the galaxy in the 4.5 /mi image, 
the alignment is better than 1/2 of the resampled IRAC pixel, i.e. < 0.2". Cyg A was only slightly 
saturated in the 8.0 /mi image but inspection of field sources appearing in both IRAC channels 
show the registration to be similarly good. 

We performed aperture photometry using 4.88" (12 resampled pixels) diameter circles cen- 
tered on the hotspots and identical adjacent apertures for background determination (Figure 2). 
The aperture correct ion for this aperture is 22.1% and 57.1% at 4.5 and 8.0 /mi, respectively 
(jLaine fc Reach! 120051 ) . Additional extinction correction s of 2.67% (4.5 iim) and 2.31% (8.0 /xm) 
were estimated from the B-band extinction of 1.644 mag (ISchlegel et alJll99Sf) extrapola ted to our 
bands assuming A^^/Ab = 0.0174 and ^.O/huMb = 0.0151 (jRieke fc Lebofskvill985l ). 



2.2. Radio Data 



There is a wealth of information on the Cyg A hotspots in the literature. We focused on 
compiling radio interferometry data with sufficient resolution and sensitivity to detect and separate 
the differing contributions of the primary (i.e., weaker and compact) hotspots B and E from the 
secondary (i.e., brighter and extended) hotspots A and D. IWright fc Birkinshawl (J200J) reported 
such a study utilizing the Very Large Array (5 and 15 GHz) and the BIMA Array (87 and 230 GHz) 
data and we have used these measurements. These authors used 10% calibration uncertainties 
for the VLA data and 12% for the BIMA ones. The BI MA 230 GHz images supersede previous 
non-interferometric mm/submm-wave measurements (e.g.. iRobson et al.lll998l ). In particular, the 
previous data were of lower resolution and did not resolve the primary from the secondary hotspots 
so did not suit our purposes. Add itionally, the integrated fluxes of the SE and NW hotspots in the 
SCUBA data (IRobson et al.lll998l ) are systematically a factor of ~two larger than the equivalent 
BIMA fluxes. We do not have an explanation for the discrepancy, but simply comment that the 
BIMA data are more consistent with the lower frequency (5 — 87 GHz) extrapolation (Figures 3-6). 



At frequencies below 1GHz, it was noted by ICarilli et al.l (|199ll ) that the spectra of the sec- 
ondary hotspots A and D flatten. The reason for such could be intrin sic flattening or cu t-off's in the 



hotspots' electron energy distribution at low energies (as favored bv lCarilli et al.lll99ll ). or an ab 



sorption effect, either synchrotron self- absorption or free- free ab sorption by extern al (foreground) 
thermal medium. The most recent high-resolution analysis by lLazio et al.l (J2006J) confirmed the 
presence of low-frequency curvature, indicating in addition that some kind of absorption (being 
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pronounced below 100 MHz) seems to be present in hotspot A, but not necessarily in hotspot D. 
For this reason, we restrict our discussion to data at frequencies > 100 MHz for which absorption 
effects are expected to be minor (if present at all), and the interferometers' beams are small enough 
to extract fluxes from the same hotspot volumes as specified in the previous section. In particular, 
we analyzed the 327 MHz map ( VLA + Pie Town Link Connected Interferometer) with a 2.5" beam 
kindly provided by T. J. W. Lazio, and measured the fluxes for hotspots A, B and D as listed in 
Table 1. Si milarly, we utilized a 1.34 GHz VLA image (provided by C. L. Carilli, and discussed 
previously in ICarilli et al.1 Il99ll ) . 



2.3. Previous Optical and Infrared Observations 



The optical and infrared observations of Cyg A are hampered by the fact that this source is at 
low Galactic latitude (6 = 5°). In addition, bright stars and other optical foreground/background 
objects are present in the close vicinities of the Cyg A hotspots, making it difficult to precisely 



determ ine the intensity of the relatively weak hotspot emission. Because of this. lMeisenheimer et al 



(|1997l ). using the United Kingdom Infrared Telescope, were able to provide only upper limits (see 
Table 1), for example < 50 ^Jy and < 5 //Jy for hotspots A and D, respective l y, at 0.452 um. 
Similarly, studies with the Infrared Space Observatory reported by ICarilli et al.l (|1999l ) resulted 
in upper limits for the hotspots' emission at 12 fim, indicating however some attenuation of the 
hotspots' synchrotron fluxes at infrared wavelengths when compared with the extrapolation of 
their radio continua. Here we adopt conservatively the 3a upper limits from ISO observations, 
namely 285 //Jy and 255 /iJy for hotspots A and D, respectively (see Table 1). We also note that 



these ISO limits are above the IRAC extrapolations for all the discussed features. iNilsson et al 



(|1997l ) performed careful analysis of the archival Hubble Space Telescope data for Cyg A, and found 
non-thermal optical emission associated with hotspots B and D. The work by lNilsson et al.1 (|1997l ) 
remained almost completely unnoticed, but here we include their measurements (converted from 
magnitudes to the flux units as given in Table 1), arguing that they match the broad-band spectra of 
both hotspots. In particular, we show that these optical fluxes are consistent with the low- frequency 
segment of the SSC component. 



2.4. X-ray Data 



The X-ray fluxes in this paper are taken from lWright &: Birk inshaw (2004), who re-analyzed 



the archival Chandra data for Cyg A reported previously by 



Wilson et al 



(|2000h . In particular, 



we adopt the results of their modeling involving an additional intrinsic absorbing column density 
for all the X-ray detected hotspots A, B, and D, being in a range = (3.6 — 4) x 10 21 cm~ 2 , in 
addition to the canonical Galactic absorbing column toward Cyg A source Ah = 3.06 x 10 21 cm -2 
(see Table 1). The additional X-ray absorption leads to slightly larger intrinsic X-ray fluxes and 
slightly steeper 0.5 — 6keV intrinsic continua for the discussed features, than obtained previously. 
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3. Modeling the Broad-Band Emission 

Figures 3 and 4 show the broad-band spectra and spectral energy distributions, respectively, 
for all the Cyg A hotspots. All of these spectra at low- frequencies, v < 1 GHz, can be fitted by 
flat power-laws with spectral indices a < 0.5 (where S u oc v~ a ), except for the weakest hotspot 
E for which the 327 MHz flux could not be determined. At higher frequencies, v > 1 GHz, the 
spectra steepen up to a > 1, and cut-off or break again around > 10 fim as indicated by the Spitzer 
detections (and also ISO upper limits). These are then the synchrotron components of the Cyg A 
hotspots. The optical < 1 /j,m emission detected from hotspots D and B is stronger than expected 
from extrapolation of the synchrotron continua to higher frequencies in the case of hotspot D, but 
not necessarily in the case of hotspot B. As mentioned above, this emission is however consistent (in 
both cases) with the low-frequency segment of the SSC component detected at keV photon energies. 
Moreover, the optical-to- X-ray power-law slopes ao-x ^ 0.6 are in reasonable agreement with the 
X-ray spectral indices found in the previous studies. We note that there is no other hotspot for 
which SSC emission has been detecte d at optical frequencies, with the possible exception of the 



radio galaxy 3C 196 (IHar dcastld 1 200 ll ) . for which however no X-ray or infrared data are available. 



In general, hotspots A and B on the jet side are much brighter in the infrared than hotspots D 
and E on the counterjet side. This is most clearly manifested for the secondaries A and D, which 
are comparably bright in radio (especially at low frequencies) and in X-rays, but which differ more 
than an order of magnitude in flux at 4 fim. The infrared continuum of hotspot A seems to be 
flatter than that of hotspot D, while at (1 — 100) GHz frequencies the opposite behavior can be 
noted. Possible reasons for these are discussed later on. Meanwhile, we concentrate on modeling 
of the multiwavelength emission of the bright hotspots D and A. 



3.1. Hotspot D and A: SSC Emission and Energy Equipartition 

Comparing the equipartition value of the magnetic field intensity deduced from the observed 
synchrotron emission of the non-thermal magnetized plasma, B eq , and the value suggested by 
the analysis of the accompanying SSC radiation sampled by a single (e.g., X-ray) flux density, 
-Bi c , is never a straightforward procedure. One of the problems here regards evaluation of the 
minimum energy condition, for which the crucial issue is the detailed shape of the synchrotron 
continuum, being poorly constrained by observations in some segments. Thus, more or less arbitrary 
assumptions have to be invoked in this respect, which may influence the resulting value of B eCL . 
The other problem is connected with the very nature of the SSC emission, namely with the fact 
that in inverse-Comptonization of the broad-band target photon population by the broad-band 
electron population there are particles and target photons with different energies which contribute 
to the monochromatic SSC flux at some particular inverse-Compton photon energy. As a result, 
single-flux-estimated B\ c is in fact a function of the partly assumed synchrotron spectral shape (or, 
equivalently, of the electron energy distribution). However, relatively well-sampled synchrotron 
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continua of hotspots D and A, together with the SSC component detected at two (optical and 
X-ray) frequencies in the former case, allow us to minimize the effects of the assumptions. 



We model the synchrotron component of hotspot D as 



-a 1 



,—Ct2 



exp(-i//i/ n 



for 
for 



n < V < U c 



(1) 



and Sl jn oc z// 3 for v < f mm . Here ct\ = 0.21, a.2 = 1.1, f cr = 3 GHz, and f max = 0.9 x 10 13 Hz 
are the spectral indices and critical (break) and maximum frequencies implied by the radio-to- 
infrared data as given in Table 1, while z/ m j n = (3 e B/Air m e c) 7^ in « 4.2 x 10 6 (B/G) 7 min Hz is the 
synchrotron frequency corresponding to the magnetic field intensity B (measured in Gauss units) 
and the minimum electron Lorentz factor 7 m i n . We treat B and 7 m i n as free parameters. For such 



m 



Carilli et al 



1991 



Meisenheimer et al. 



1997 




Wilson et al. 


2000; 



Kino fe Takaharal 12004 ). we 



construct the implied electron energy spectrum, n e {j) for 7 m i n < 7 < 7 max , with the normalization 
K e and the appropriate critical energies 7 cr and 7 max depending on B (see § 4 and equation 3). 
Next, we compute the SSC emission using the useful approximation formula 



Sf, 



<7T B.V 



1/2 



3i//4 7] 
3y/4 7] 



2 

min 




(2) 



(see, e.g., IChiang k, Dermerlll999l ). where R = 0.8 kpc is the effective radius of the hotspot, i.e. 
radius of a sphere with volume equivalent to the cylindrical structure measured directly from the 
radio maps. All the model parameters and model results discussed in this section are summarized 
in Table 2. 

The evaluated SSC emission at 1 keV photon energies is consistent with the observed X-ray flux 
of hotspot D for B lc = 270 /uG and any minimum electron Lorentz factor 7 m i n lower than or equal 
to 300. We treat the latter value as the upper limit for 7rn i p, sin ce Cyg A hotspots were detected 
at < 100 MHz frequencies (jCarilli et al.Ml99ll : lLazio et al.1 120 06). while synchrotron radiation of 
7 = 300 electrons in B = 270 fiG magnetic field corresponds almost exactly to v = 100 MHz. 
The resulting broad-band synchrotron and SSC emission of hotspot D is shown in Figure 5 (solid 
lines). We note, that the obtained magn etic field intensity is in the higher end of (although still 
consistent with) the previous es timates ([Harris et al.l 11994 ; IWilson et al.l [200CJ; iKino &; Takahara 
20041 ; IWright &: Birkinshawil2004l ) . This agreement — despite the slightly different electron spectral 
shape considered by different authors — is due to the fact that for the hotspot D parameters, 
electrons with energies 7 ~ 10 4 , i.e. the ones emitting relatively well studied ~ 100 GHz synchrotron 
emission, are primarily responsible for production of the X-ray SSC flux. It is worth mentioning 
that our mode ling yields also a value of Qjg> c y V = 0.93, which is comparable to the X-ray spectral 
index given by lWright &: Birkinshawi ([200J), namely ax = 0.8 ± 0.11. 



Interestingly, with B lc = 270 fiG the expected optical-to-X-ray power-law slope of the SSC 
emission changes slightly from ao-x = 0.65 for 7 m i n = 1 down to «o-x = 0.6 for 7 m j n = 300. These 
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values are all consistent with the observed one ao-x ~ 0.6, and the flux found by iNilsson et al 



(jl997l . as a residual flux after subtraction of a bright star in the field) is overproduced by no more 
than a factor of 1.8. This is in fact a crucial result, since the optical SSC emission is produced 
mostly by the electrons with energies 7 < 10 3 , and thus probes the low-energy segment of the 
electron distribution. Indeed, if one assumes a standard i/~ 5 synchrotron spectrum below the 
break frequency v cr instead of the flat cx i/~ 0,21 continuum considered above (arguing that the 
observed flux at low radio frequencies may be significantly modified by absorption effects), the 
implied magnetic field intensity needed for producing the X-ray SSC emission would be almost 
unchanged (namely B 1C = 290 fJ.G), but the optical emission would be overproduced by a factor of 
> 3 for 7 m j n = 1—300 (see dotted lines in Figure 5). Thus, the flat synchrotron continuum of hotspot 
D below f cr seems to be intrinsic to the source, and not resulting exclusively from absorption effects. 
Note in addition, that the choice ot\ = 0.5 would lead to overproduction of the 4.5 /jm emission 
as detected by Spitzer. Of course, inevitable synchrotron self-absorption occurring at < 100 MHz 
frequencies is expected to decrease the predicted SSC flux at low (< 10 15 Hz) photon energies, and 
such a decrease should be stronger for steeper synchrotron spectra at low frequencies. However, 
even in the case of a\ = 0.5 this effect would not remove the aforementioned overproduction of the 
infrared and optical fluxes of hotspot D. 

Finally, we also check if possible steepening of the synchrotron continuum in the 'unconstrained' 
frequency range > 230 GHz could somehow change the conclusions presented above. In particular, 
we repeat the SSC modeling of hotspot D as discussed previously but with an additional spectral 
break by Aa = 0.5 at v\> T = 0.5 x 10 12 Hz. In § 4.2 we argue that around this frequency — for the 
particular parameters of hotspot D considered — one should indeed expect radiative cooling effects 
leading to steepening of the synchrotron spectrum. The analysis presented in Figure 5 (dashed 
lines) indicates however that the introduced additional spectral break affects only the unobserved 
high-energy (^$> 1 keV) part of the expected SSC emission, but not the derived value of B lc nor the 
low-energy tail of the SSC component. 

We follow up the SSC analysis for the case of hotspot A, assuming in a first approach a single- 
broken spectral shape of its synchrotron continuum. From the collected radio-to-infrared data we 
infer the input parameters a± = 0.28, a<i = 1.2, v CT = 2.6 GHz, and R = 1.1 kpc (see Table 2). 
Regarding the high-energy cut-off, we find that the value i^ max = 3.25 x 10 13 Hz leads to a rough 
agreement with the detected 4 — 8 /jm flux, although the implied synchrotron spectrum is then 
inconsistent with the relatively flat infrared slope constrained by Spitzer, and it also violates the 3a 
upper limits by ISO. In other words, the infrared spectrum of hotspots A cannot be really modelled 
in terms of high-energy synchrotron cut-off, contrary to the case of hotspot D. Slightly better fit 
is however obtained by introducing an additional break (as before by about Aa = 0.5) around 
z^br = 1-2 x 10 12 Hz (see § 4.2), followed by an exponential cut-off at ^ max > 10 14 Hz. Note, that in 
such a case the maximum synchrotron frequency cannot be constrained with the available dataset. 
With the two aforementioned descriptions of the synchrotron spectrum, we find that in both cases 
the expected SSC radiation matches the observed X-ray flux of hotspot A for B 1C = 170 fiG and a 
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broad range for 7 m i n = 1 — 300 (see solid and dashed lines in Figure 6). Clearly, the two different 
cases considered are energetically equivalent. The evaluated power- law slope Qjg^ y V = 1 is now 



slightly larger than the appropriate X-ray spectral index provided bv lWright Birkinshawl (|2004l ). 
ax = 0.77 ±0.13. Again, the assumed i^ -0 5 synchrotron spectrum below u cr leads to slightly larger 
value of B lc = 175 /uG, but this time to the optical/infrared fluxes well below the available upper 
limits/detections, and so cannot be formally excluded (see dotted line in Figure 6). 

Figure 7 shows the ratio of energy densities in the radiating electrons and magnetic field for 
hotspots D and A inferred from SSC modeling, U e /UB, as a function of the unknown minimum 
electron Lorentz factor 7 m i n . One can see that in the case of hotspot D (lower solid line) electrons 
dominate energetically over the hotspot's magnetic field by only a small factor of 4 — 3 for the 
considered range 7 m i n = 1—300. In the case of hotspot A (upper solid line) this effect is bit stronger, 
since U e /U-Q ~ 8 — 6 (see also Table 2). Moreover, the oc z/~ ' 5 low-frequency synchrotron spectra 
considered previously for both hotspots as possible (though not likely) case would result in a larger 
deviation from the minimum energy condition unless 7 m i n > 100 (lower and upper dotted lines for D 
and A, respectively). Thus, we conclude that the second ary hotspots in Cyg A seem to be particle 



rather than magnetic field dom inated (see in this context IWilson et al 



2000 



2004 



Wright &: Birkinshaw 



Kino Takaharal [2004J) . However from one perspective the derived energy density ratios 



are surprisingly close to equality. For example, if some fraction of the o bserved X-ray emission 
were to come from some other emission mechanism (e.g., synchrotron; see lHardcastle et al.l [200J ; 
Balucinska-Church et al.ll2005l ) then the calculated value of B 1C would become a lower limit on B. 
This would mean fewer electrons to produce the observed radio emission and thus the value of 
U e /U-Q would be reduced. Other conditions such as a filling factor < 1 would also lead to stronger 
fields and fewer electrons. 



3.2. Hotspots/Counter-Hotspots Infrared Asymmetry 



The fact that the hotspots A and B on the jet side are much brighter in infrared than hotspots 
D and E on the counterjet side, while being roughly comparable at other wavelengths (in particular 
secondary hotspots A and D), raises the question on the nature of the observed asymmetry. There 
may be several possible reasons for this. One is that due to the intermittent nature of the jet 
terminal regions and light-travel effects, hotspots on the jet side are observed at a different phase 
of their evolution than the an alogous features on the counterjet side. For the jet viewing angle ~ 
60° — 70° in Cyg A source (see lCarilli Barthellll996l ). the jet/hotspot intermittency on timescales 
< 10 5 yrs would be then req uired to observe any differences in the hotspots' appearance. In fact, 
Wright Birkinshawl (120041 ) analyzing the radio data noted that 'less spectrally aged material is 
radiating in hotspot D than in A', consistent with the idea that the observed secondary hotspot on 
the jet side is older than its counterpart (but see § 4.2). Yet it is not clear why the more evolved 
feature should be characterized by higher maximum synchrotron frequency, as suggested by the 
Spitzer observations reported here. 
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Another possibility is that the time-travel effects/jet intermittency does not play any role, 
but that instead the electron acceleration process at terminal shocks acts a bit differently on 
opposite sides of the nucleus. The reason for this could be for example a slightly different con- 
figuration of the magnetic filed with respect to the shock front, assuming that the upstream (jet) 
plasma is moving with at lea st mildly relativistic velocities. We note, that numerical simulations by 



Niemiec &: Ostrowskil (120041 ) clearly indicated that the particle spectra (in particular their slopes 
and maximum energies) resulting from the diffusive (1st order Fermi) acceleration process are in- 
deed very sensitive to the assumed magnetic and kinematic parameters of relativistic shocks. This 
possibility could be studied by detailed analysis of the Cyg A hotspot regions at radio frequencies, 
namely by investigation of the hotspots' polarization properties. We note in this context, that our 
modeling indicates some differences in both magnetic field intensity and high energy synchrotron 
cut-off between hotspots A and D. This suggests in turn magnetic-related differences in the particle 
acceleration processes at the highest (7 ~ 10 5 and above) electron energies. 

An interesting issue to note is that the obtained differences in the magnetic field intensity B 1C 
and in the equipartition ratio U e /U-Q for hotspots A and D do not imply different total kinetic 
luminosities of the radiating plasma within these two features. Indeed, the total hotspots' energies 
can be estimated as E tot = (U e + Ub) x V, where V = |vri? 3 . The SSC modeling allows us to 
find E to t ~ 1.7 x 10 57 erg for hotspot A, and E to t ~ 0.9 x 10 57 erg for its counterpart D (assuming 
7min = 1)- Meanwhile, the hotspots' dynamical timescale can be evaluated simply as idyn = 2R/f3c, 
where (3c is the bulk velocity of the emitting matter. Taking (3 = 1/3 as appropriate for the plasma 
downstream of a strong relativistic shock front, one can find id yn ~ 2.2 x 10 4 yrs and id yn ~ 
1.6 x 10 4 yrs for A and D, respectively. Thus, the implied total kinetic luminosity L^ m = E tot /td y n 
is very similar for the two analyzed features, being roughly L^m ~ 2 x 10 45 erg s _1 . We also note 
that twice this luminosity — i.e., total power carried by the radiating electrons and magnetic field 
in both hotspots — is only slightly lower than total kinetic power inferred from expansion of the 



Cyg A cavity in the surrounding cluster gas, ~ 10 46 erg s" 1 (jWilson et al.1 12009 1. 



4. Constraining The Electron Energy Distribution 

Electron spectra of hotspots D and A within the energy range 300 < 7 < 30000 (not af- 
fected by absorption or cooling effects), as found from the SSC modeling presented above, can be 
approximated by a broken power-law 



/ \ r , j 7 p for 7 < 7 cr 

n c 7 =K c x ■> , 3 

7cf 7 v for 7 > 7 cr 



with 7 cr 2 x 10 3 , pi = 2 a\ + 1 rj 1.4 — 1.6, and P2 = 2 02 + 1 ~ 3.2 — 3.4. At higher energies the 
electron distributions — possibly modified by the radiative loss effect — seem to continue up to at 
least 7 ~ 10 5 . These spectra are shown in Figure 8 (solid and dashed lines correspond to hotspots 
D and A, respectively; see also Table 2). We emphasize that 



-11- 



(i) . for both hotspots the critical break energy corresponds almost exactly to the mass ratio 

between protons and electrons, 7 cr ~ m p /m e , which gives a natural scale unit in the case of 
a shock front formed in an electron-proton plasma; 

(ii) . for both hotspots the low-energy segment of the electron distribution continues down to at 

least 7 m i n ~ 0.1 m p /m e electron energies, with flat power-law spectrum of slope p\ « 1.5 
possibly slightly different in the two analyzed cases; 

(hi), for both hotspots the high-energy part of the electron continuum is very steep, with spectral 
indices p2 > 3 and maximum energies 7max ^ 50 m p /m e slightly different for the two features; 

(iv). in both hotspots 'standard' electron spectrum n c (^) oc 7~ 2 expected from the diffusive (1st 
order Fermi) shock acceleration in the non-relativistic test-particle limit is not observed. 



The latter finding regarding the model shock-type spectrum oc 7 -2 is an interesting issue. It should 
not be surprising though, keeping in mind the fact that the terminal shocks of powerful jets, like 
those of the Cyg A radio galaxy analyzed here, are almost for certain (mildly) relativistic, with 
an oblique magnetic field configuration. As such, they are indeed not expected to be the sites 
of the diffusive accel eration process known from the non-relativistic test-particle mo dels (see the 
relevant discussion in lBegelman &: Kirkl|l990l ; lNiemiec et al.ll2006l ; iLemoine et al.ll2006l ). The break 
frequency corresponding to the m p /m e mass ratio preceded by a flat-spectrum power-law is also an 
interesting but again not a very unique (see lLeahy et al.lll989l ) nor completely unexpected finding. 
As discussed below, it may imply several important constraints on extragalactic jets in general. 



4.1. Low-Energy Spectrum 



The question of accelerating the electrons from thermal (non- or mildly-relativistic) to ultra- 
relativistic energies at relativistic shocks is widely debated. That is because the electrons have 
to be already ultrarelativistic to undergo diffusive lst-order acceleration, since otherwise their gy- 
roradii are smaller than the thickness of the velocity transition region, and thus the low-energy 
particles c an not cross the shock fr ont freely. This issue, known as the 'injection problem' (e.g., 
Bellll978al lbl:l Eilek &: HughesHl990l ). is of particular relevance for the collisionless plasma which is 
dynamically dominated by protons, because the shock thickness is then expected to be of the order 
of the dominant protons' gyroradius. For the upstream matter being cold in its rest frame, the 
required 'injection' electron Lorentz factors would be roughly ~ m p /m c . It was postulated, that 
some efficient but unspecified energy/temperature coupling acting between protons and electrons 
(which cannot be collisional in nature) should serve as the pre-acceleration mechanism. Yet the 
low-energy (7 < m p /m e ) electron spectrum resulting from such a process remai ned unknown. The 
Cyg A hotspots are, o bviously, subject to this problem (see the discussion in ICarilli et al.lll99ll ; 
Kino k Takaharall2004l l. 
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In their pioneering work, iHoshino et al.l (|1992l ) presented a possible solution to the injection 
problem for the case of a cold electron-proton plasma upstream of the shock. In particular, by using 
ID particle-in-cell simulations, they showed that within the velocity transition region positrons can 
absorb electromagnetic waves emitted at high harmonics of the cyclotron frequency by the cold 
protons reflected from the shock front, and thus be accelerated to suprathermal energies. For 
efficient acceleration it was required to assume at least comparable energy densities in protons and 
in electron-positron pairs. T his result conce rned strictly positrons and not electrons, because for 
the parameters considered by lHoshino et al. polarization of the pro ton cyclotron waves was almost 
exclusively left-handed. Recently, however, lAmato Aronsl (|2006l ) — who applied more r e alistic 
parameters in their numerical studies when compared with simulations by IHoshino et al.l ([1992) 
- showed that the same process may as well apply to the electron^. The resulting electrons 
spectra are consistent with a power-law between energies 7 ~ T, where T is the bulk Lorentz 
factor of the upstream medium, and 7 ~ T (m p /m e ), possibly reduced by a factor of a few due 
to thermal dispersion in the upstream proton momenta. The spectral index of this distribution, 
hardly constrained by the simulations, seems to depend on the plasma content (i.e., on the ratio of 
proton number density to the electron number density ) , and can be relatively flat. We em phasize 
that although simulations presented by IHoshino et al.1 (|1992l ) and lAmato Aronsl (|2006l ) involve 
large bulk Lorentz factors of the upstream plasma, analogous processes are expected to take place 
also at mildly relativistic shocks considered in this paper (T ~ Tj ~ few), so the model results can 
be qualitatively applied. 



Our analysis indicates that a flat power-law electron spectrum oc 7~ L5 in Cyg A hotspots ex- 
tends from low (unfortunately unconstrained) electron energies up to almost exactly j CT ~ m p /m e . 
This is in fact in a good agreement with t he simple (ID, ultrarelativistic) resonant acceleration 
m odel discussed by lAmato Sz Aronsl (|2006l ). If such an association is correct, and if the results 
of lAmato k, Aronsl as applicable, then it would automatically imply that (i) the jets in powerful 
radio sources like Cyg A are made of electron-proton rather than electron-positron plasma, (ii) 
a significant fraction of jet kinetic power is carried by cold protons, and (in) number density of 
protons within the jets is most likely lower than that of the electron-positron pairs. We note, that 
th e above conclusions regard ing the jet content would be then consistent with the ones presented 
by ISikora &: Madeiskil (|2000l ). who analyzed the broad-band spectra and vari ability of radio-loud 
quas ars, believed to be beamed counterparts of FR II radio galaxies (see also iGhisellini h Celotti 



20021 . who presented other evidence for the dynamical role of protons in quasar jets) 



2 For the parameters used in the ID particle-in-cell simulations reported by Hoshino et al. 1 19921 ) — small proton- 
to-electron mass ratio m p /m c ~ 10 and similar energy densities in leptons and hadrons — the number of protons was 
comparable to the number of pairs, and thus polarization of the protons' cyclotron emission was almost completely 
left-handed. Amato fe Aronsl 1 20od ) considered much larger proton-to-electron mass ratio m p /m c ~ 100 which, for 
comparable energy densities in the two plasma species, corresponds to a much smaller number of protons than pairs. 
Under such conditions, proton cyclotron emission is expected to be a mixture of left- and right-handed modes, which 
then can be absorbed by both electrons and positrons. 
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4.2. High-Energy Spectrum 



It was stated in previous studies (jCarilli et al.lll99ll ; iMeisenheimer et al.lll997l ). that radiative 
cooling effects are of importance for the electrons emitting > 1— 10 GHz synchrotron photons within 
the Cyg A hotspots A and D, and hence that the spectrum of 'freshly' accelerated electrons in the 
appropriate energy range is flatter by some Ap than the observed (cooled) one. In a framework 
of the 'continuous injection' model (widely applied to hotspots in general in spite of the obvious 
fact that the CI model rests on the premise that particles accumulate in the emitting volume 
whereas hot spots are emitting volumes characterized by both a hefty s ource function and a strong 
outflow; see iHeavens fc Meisenheimerl 119871 : IMeisenheimer et al.lll989l ). such a corrected electron 
spectral index would be then p2 = 2.2 — 2.4 instead of the observed one p2 = 3.2 — 3.4, assuming 
a homogeneous distribution of the magnetic field within the emitting volume leading to Aa = 0.5. 
If this is the case, then fine-tuning of the model parameters would be required in order to equalize 
the observed critical energy j CT ~ m p /m e , argued above to be intrinsic to the electron spectrum, 
with the break energy resulting from cooling effects, denoted here as 7b r - This, as discussed below, 
is however not expected. 

The synchrotron break frequency corresponding to the transition between slow- and fast- 
cooling regime, = (3 e B/Att m e c) 7^, can be found by equalizing the dynamical timescale in- 
troduced previously, idyn = 2R/f3c, with the energy-dependent radiative cooling timescale i ra d(7)- 
The latter one can be found as t ra d = 7/|7|syn> where the (dominant) synchrotron cooling rate is 
simply |7|syn = cot B 2 j 2 / (6ir m e c 2 ) . Hence, the condition td yn = *rad(7br) gives 



^br 



277r m e c 3 ef3 2 
44 B 3 R 2 



6 x 10 13 (3 2 Bzl R^ c 



Hz 



(4) 



where 5_ 4 = S/10" 4 G, and R kpc = R/l kpc fsee lBrunetti et al.ll2003l : ICheung et ajj2005j ). Taking 
hotspots' parameters as discussed in previous sections, one obtains = 0.5 x 10 12 Hz for hotspot 
D and z/b r = 1-2 x 10 12 Hz for hotspot A. This is indeed much higher than the observed critical 
frequency v CT ps 3 x 10 9 Hz. Hence we conclude that the observed steep electron spectrum at 
7 > m p /m e results directly f rom the acceleration pro cess, and is not affected by the radiative 
cooling effects. We note that iKino Takaharal ((2004J) evaluated the cooling break frequency a s 
hig h as discussed above. The d iscrepancy with the previous works, like that bv lCarilli et al.l (|199ll ) 
and IMeisenheimer et al.l (|1997l ) is mainly due to a much lower outflow velocity assumed by these 
authors for the emitting region downstream of the terminal shock front. This resulted in a larger 
(by more than an order of magnitude) dynamical timescale for the hotspot, and thus a lower cooling 
break frequency. From equation 4 it follows directly that the outflow velocity should be reduced as 
much as (5 < 0.03 in order to obtain v^ r < 10 GHz. Even keeping in mind the expected very complex 
structure of the outflowing plasma in the discussed region, such low downstream velocities seem to 
be not realistic. Note, that typical advance velocities of hotspot features in FR II radio galaxies 
which have to be lower than (3 discussed here — are typically /3 ac i v 0.01 (e.g.. iMachalski et al 



20071) 
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If we consider exclusively the case where diffusive shock acceleration is responsible for the 
formation of the high energy electron spectra observed in the Cyg A hotspots, we have to accept 
a significant departure of the electron spectral index from the value P2 = 2 corresponding to the 
non-relativistic shock, but also from the often claimed 'universal' value P2 = 2.2 corresponding to 
the (ultra) relativistic shock. In fact, recent mod e ling of the first-orde r Fermi acceleration process 
at relativistic shocks (INiemiec &; Ostrowskil I2004J : iNiemiec et al.ll2006l ) indicates the possibility for 
formation of extended high energy 'power-law + cut-off' tails above the shock-related F m p c 2 energy 
scale. The theory is however not able to provide robust evaluations of the involved spectral indices 
and maximum energies yet, and thus the analysis presented here could be considered as valuable 
empirical measurements of these parameters for mildly relativistic shocks. We emphasize that 
the observed differences in the power-law slopes and cut-off frequencies between the two hotspots 
may then indicate sensitivity of the particle spectra resulting from the 1st order Fermi acceleration 
process on (even relatively small) differences in the magnetic field configuration and intensity within 
the shock front. We note in this context, that hotspot A shows a larger deviation from the minimum 
power condition than hotspot D. 

The second possibility we consider is that of the shock responsible for accelerating the elec- 
trons only below ~ T m p c 2 energies (as discussed in § 4.1), and of a distributed acceleration in the 
turbulent hotspot volume responsible for accelerating electrons above that scale. The acceleration 
processes involved can be the second-order Fermi acceleration by relativistic hotspot turbulence 
downstream of the jet terminal shock, or particle acceleration accompanying magnetic field recon- 
nection events. Both these processes were hardly studied for the particular hotspots' conditions, 
but the appropriate rough time scale evaluations easily allow for electron energies as high as the 
7max values required here. A s illust r ative examples o f recen t theoretical progress in this field one 



can compare^ , , >,..„,„ v . , v ..... -> ..*..„... 

.laroschek et al ( 20041 ): Lyubarsky ( 2005 ): Zenitani Hoshino ( 2005 ) for the reconnection process 



Petrosian fc Liul ( 2004 



Cho & Lazarian (2006) for the second-order acceleration, and 



Also in this case the electron spectral slope derived from the observations could serve as an impor- 
tant input for better understanding of the aforementioned acceleration processes. 



5. Conclusions 

We have presented the results of Spitzer Space Telescope observations of the famous radio 
galaxy Cygnus A. We have detected both hotspots A and B on the west (jetted) side at 4.5 fim 
and 8 |im. Because confusing emission precludes accurate estimates of the background level, our 
measurements of D and E on the eastern side are of low statistical significance. In all cases the 
observed infrared emission is consistent with the high-energy tail of the hotspots' synchrotron radi- 
ation. When combined with the other data collected from the literature, our observations allow for 
detailed modeling of the broad-band (radio-to- X-ray) emission for the brightest hotspots A and D, 
uniquely sampled at different wavelengths. We confirm that the X-ray emission detected previously 
by ROSAT and Chandra satellites is consistent with the synchrotron self-Compton radiation for 
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the magnetic field intensity B ~ 170 fiG in hotspot A and B ~ 270 //G in hotspot D. However, we 
also find that the energy density of the emitting electrons is larger than the energy density of the 
hotspots' magnetic field by a factor of a few, and that hotspot A is more particle dominated than 
hotspot D. Finally, we note that the hotspots located on the jet side (A and B) are significantly 
brighter in the infrared than their counterparts located on the other side of the active center (D and 
E). We suggest that such an asymmetry is most likely due to different configurations and different 
relative strengths of the magnetic field within the terminal shocks on the jet and counterjet sides. 
These, in turn may result either from temporal evolution of the hotspots (coupled with the fact 
we are observing them at different proper ages), from differences in the cluster medium on the two 
sides, or from differences in the jet energy or field configuration on timescales of 10 5 yrs. 

The performed analysis allows us to constrain relatively precisely the energy spectra of ultra- 
relativistic electrons producing the broad-band emission in hotspots A and D. We find that in both 
cases it is consistent with a flat power-law n e (j) oc 7 _Pl with pi = 1.4— 1.6 for the electron energies 
7 < m p /m e , followed by a steep power-law n e (7) oc r y~ P2 with p2 = 3.2 — 3.4 for 7 > m p /m e . We 
argue that such a shape of the electron continuum reflects most likely two different regimes of the 
electron acceleration process taking place at mildly relativistic shock, rather than resulting from 
the radiative cooling and/or absorption effects. In this picture the protons' inertia defines a critical 
energy for the hotspot electrons above which Fermi-type acceleration processes may play a major 
role, but below which the operating acceleration mechanism has to be of a different type. We 
suggest that the latter mechanism is connected with collisionless processes acting within the shock 
transition layer, most likely involving; resonant absorption of the shocke d ions' cyclotron waves by 
the electron-positron pairs (|Hoshino et al.1 Il992l : lAmato &; Aronsl 120061 ) . If this interpretation is 
correct, it would automatically imply a dynamical role of (cold) protons within the jet, consis- 
tent with analyses of jet energetics in quasar sources. It is important to note that the hotspots 
in other powerf ul radio sources analogous to Cyg A reveal similar spectral breaks around GHz 
frequencies (see iLeahv et al.l Il989l ) . Obviously, insufficient understading of the relativistic shock 
structure, as well as limitations of the available numerical methods in modeling microphysics of rel- 
ativistic magnetized plasma, leave some room for the other processes which can be responsible for 
(pre) acceleration of low-energy electrons. Let us mention in this context the in teraction of particles 
with electrostatic waves excited at the front of quasi-perpendicular shocks (jShimada Hoshino 
2000 



Hoshino Shimada 



stream shock instabilities ( 



2002 ), or with quasi-stationary electric fie ld associated with the two 



Hededal et al.l 12004 : iNishikawa et al.l 120061 ) . We note that any type of 



turbulence inevitably generated in e _ p + shock by protons will lead to acceleration of electrons up 
to the downstream proton energies, but only a few simple models regarding this type of processes 
were studied till now. 

The high-energy part of the electron distribution above m p /m e energies may result from the 
Fermi-type processes, although the steep slope of the electron spectrum at these energies deviates 
significantly from the 'universal' spectral indices often claimed in the literature for diffusive shock 
acceleration. We argue however that such a disagreement is in fact only superficial, since the most 
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recent theoretical studies on the particle acceleration taking place at relativistic and oblique shocks 
clearly indicate an absence of any universal spectrum, but instead show a variety of particle spectral 
shapes. This is supported further by the observed differences in the electron spectral shape between 
hotspots A and D at the highest electron energies, > 50 m p /m e , which indicates that the 1st order 
Fermi acceleration process — if indeed responsible for the formation of the high-energy electron 
tail — is very sensitive to (even relatively) small differences in the magnetic field structure at the 
shock front. We also note that our conclusions are consistent with the general finding that the 
electron spectrum injected from the terminal hotspots to the lobes of powerful FR II radio galaxies 
is not of a single and universal power-law form (as discussed in detail by, e.g.. iRudnick et al.lll994l ; 



Treichel et al 



2001 



Machalski et al.ll2007l ). This, often ignored result, has several consequences for 



estimating the jet lifetimes and thus energetics and duty-cycles of radio-loud active galactic nuclei. 
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Table 1. Spectral flux densities for the Cyg A hotspots. 





Hotspot A 


Hotspot B 


Hotspot D 


Hotspot E 




v [Hz] 


Su [Jy] 


Su [Jy] 


Su [Jy] 


Su [Jy] 


ref. 


0.327 E9 


173±25 


50±15 


162±20 




[1] 


1.34 E9 


116.4±5.8 


31.5±3.0 


120.6±6.0 


10±1.0 


[2] 


5.0 E9 


44.60±4.46 


9.96±1.00 


57.69±5.77 


3.90±0.39 


[3] 


15 E9 


10.80±1.08 


2.47±0.25 


17.74±1.77 


1.16±0.12 


m 


87 E9 


1.49±0.18 


0.44±0.05 


2.49±0.30 


0.19±0.02 


m 


230 E9 


0.43±0.05 


O.llrbO.Ol 


0.87±0.10 


0.10±0.01 


[3] 


2.5 E13 


< 285 E-6 


< 285 E-6 


< 255 E-6 


< 255 E-6 


[j] 


3.798 E13 


156±19 E-6 


50±12 E-6 


31±25 E-6 


15±24 E-6 


[5] 


6.655 E13 


75±9 E-6 


25.3±18 E-6 


5.5±5.5 E-6 


3.5±7 E-6 


[5] 


1.36 E14 


< 30±3 E-6 


< 35±2 E-6 


< 25 E-6 




[6] 


2.48 E14 


< 100±10 E-6 


< 52.4±2.6 E-6 






[6] 


4.29 E14 






1.63±0.3 E-6 




[7] 


4.63 E14 


< 80±8 E-6 


< 46.2±0.3 E-6 






[6] 


5.45 E14 




2.65±0.6 E-7 








6.64 E14 


< 50±5 E-6 


< 21.9±0.8 E-6 


< 5 E-6 




[6] 


2.42 E17 


31.2±4.3 E-9 


6.8±2.6 E-9 


47.9±5.9 E-9 




[8] 


«x 


0.77±0.13 


0.70±0.35 


0.80±0.11 




[8] 



(0.5-7.5 keV) (0.5-4.5 keV) (0.5-6.0 keV) 



No te. — References : [1] T .J. W. Lazio, private communicat ion ( see also lLazio et al.1 
I2006T) ; [2l lCarilli et~aH Jl99lf) ; [31 IWright fc Birkinshawl i2004) : [ 41 ICarilli et al.Nl999l. 
3a upper limits); [5] this work; [6] lMeiserAeTm^r^t^dTJl997l ); [7] iNilsson^^al* Nl997l) ; 
[8] IWright fc Birkinshawl | |2004 corrected for an additional intrinsic absorption) 
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Table 2. Model Parameters for the Cyg A hotspots 





Hotspot A 


Hotspot D 


Measured Parameters 


OL\ 


0.28 


0.21 


OL2 


1.2 


1.1 


V a [Hz] 


2.6xl0 9 


3.1xl0 9 


^max [Hz] 


>3.3xl0 13 


0.9xl0 13 


Cylindrical dimensions (L, r) 


(2.2", 0.8") 


(1.9", 0.7") 


Equivalent radius R [kpc] 


1.1 





8 


Fitted Parameters 


Tmin 


1 


300 


1 


300 


a x 


1 


1 


0.93 


0.93 


ao-x 






0.65 


0.60 


Tcr 


1.9 xlO 3 


1.9xl0 3 


1.7xl0 3 


1.7xl0 3 


Tmax 


>2.2xl0 5 


>2.2xl0 5 


0.9xl0 5 


0.9xl0 5 


K c [cm -3 ] 


1.4xl0~ 4 


1.4xl0" 4 


7.6xl0~ 5 


7.6xl0 -5 


B ic [fiG] 


170 


170 


270 


270 


U C /U B 


8 


7 


4 


3 



Note. — The spectral parameters, ai, ct2, f C r, and i/ max , are defined in equa- 
tion 1. We estimated the length L and radius r of a cylindrical volume for the 
hotspots; for convenience, we give the equivalent radii R of spheres of equal 
volume to the cylinders. 

From the SSC modeling, we infer the magnetic field B and report two spectral 
indices of the SSC spectrum which can be compared with observations (ao-x> 
and ax which is between 0.5 and 6 keV). We constrain also the parameters of 
the electron energy spectrum (7 cr , 7 max , and K e ; see equation 3), as well as the 
ratios of the energy densities U c /Ub- 
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Fig. 1. — Spitzer Space Telescope IRAC and VLA image (green; from I Per ley et al.l 1 1984 ) of 
Cygnus A. The 4.5/um sources are shown in blue, 8.0/^m emission in red and those emitting in 
both bands in pink. 
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Fig. 2. — 8^m images of Cyg A hotspots with 5 GHz contours superimposed. Regions used for the 
hotspots' and backgrounds' flux subtraction are marked by solid and dashed circles, respectively. 
Upper: Secondary hotspot A on the jet side. Middle: Primary hotspot B on the jet side. Lower: 
Secondary hotspot D (left) and primary hotspot E (right) on the counterjet side. 
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Cyg A Hotspots 




1E8 1E9 1E10 1E11 1E12 1E13 1E14 1E15 1E16 1E17 1E18 

v [Hz] 



Fig. 3. — Flux densities for Cyg A hotspots (see Table 1). Lines do not correspond to any model 
fit, but only connect data points for each hotspot. 
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Fig. 4. — Spectral energy distribution of the broad-band emission observed from Cyg A hotspots. 
Open triangles denote upper limits, open stars correspond to known foreground star in the field, 
while filled circles correspond to the detected fluxes (see Table 1). 
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Fig. 5. — Spectral energy distribution of hotspot D. Open triangles denote upper limits, while filled 
circles correspond to the detected fluxes. The parameters of the model described by the solid line 
are listed in Table 2. As discussed in § 3.1, our modeling is little changed if we input an additional 
break (u^ T = 0.5 x 10 12 Hz) in the synchrotron spectrum (dashed line). Increasing ct\ to 0.5 (from 
0.21) produces the dotted curves and overproduces the low-energy portion of the SSC component. 
All the lines assume 7 m i n = 1. 




Fig. 6. — Spectral energy distribution of hotspot A. Open triangles denote upper limits, while filled 
circles correspond to the detected fluxes. The parameters of the model described by the solid line 
are listed in Table 2. As in modeling the hotspot D (Figure 5), we show the effect of adding a break 
{iAa T = 1.2 x 10 12 Hz) in the synchrotron spectrum (dashed line), or increasing a\ to 0.5 (dotted 
line). All the lines assume 7 m i n = 1. 
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Fig. 7. — Ratio of energy densities in ultrarelativistic electrons and magnetic field for hotspots A 
(lower solid and dotted lines) and D (upper solid and dotted lines). Solid lines correspond to the 
case with low-frequency synchrotron spectra S u<Uci oc v~ ai with a\ = 0.21 (spot D) or 0.28 (spot 
A), while dotted lines correspond to ol\ = 0.5 for both features. 
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Fig. 8. — Electron energy spectrum for hotspots A (dashed line) and D (solid line), as inferred from 
the SSC modeling. Shaded regions illustrates energy ranges affected by the absorption or radiative 
losses effects. Note the 'normalizing factor' 7 2 in the abcissa scale. 



